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Abstract. We investigate dark and visible matter distribution in the Coma cluster in the 
case of the Navarro-Frenk-White (NEW) profile. A toy model where all galaxies in the cluster 
are concentrated inside a sphere of an effective radius Reff is considered. It enables to obtain 
the mean velocity dispersion as a function of Reff- We show that, within the observation 
accuracy of the NEW parameters, the calculated value of Reff can be rather close to the 
observable cutoff of the galaxy distribution . Moreover, the comparison of our toy model with 
the observable data and simulations leads to the following preferable NEW parameters for the 
Coma cluster: R 200 ~ 1-77 h~^ Mpc = 2.61 Mpc, c = 3-^4 and M 200 = 1.29/i“^ x 10^^Mq. In 
the Coma cluster the most of galaxies are concentrated inside a sphere of the effective radius 
Reff ~ 3.7 Mpc and the line-of-sight velocity dispersion is 1004kms“^. 


Contents 


1 Introduction 1 

2 Zero acceleration sphere and mean velocity dispersion 2 

3 Navarro-Frenk-White profile in the Coma cluster. Theory vs observations 

and simulations 5 

4 Conclusion 9 


1 Introduction 

According to the recent observations [1-5], our Universe is dark: dark energy and dark mat¬ 
ter contribute approximately 69% and 26% into total mass-energy balance in the Universe, 
respectively. Different independent observations also indicate that dark matter (DM) en¬ 
velops the galaxies and clusters of galaxies. Baryonic matter is only 10-15% of the total 
mass of clusters of galaxies [6]. The rest is DM. Unfortunately, the nature of DM (as well 
as dark energy) is still unclear and is a subject of numerous investigations [7]. The standard 
cosmological model, i.e. the ACDM model, assumes that DM is cold. This model is rather 
successful in explaining the structure formation in the Universe. There is a quite big num¬ 
ber of papers devoted to the investigations (analytical and numerical) of the structure and 
dynamics of galaxies and clusters of galaxies (see, e.g., the reviews [8] and [9]). For these 
investigations, the spherical density profile p{R) of DM in galaxies and clusters of galaxies is 
very useful. On the one hand, such profiles help to reveal different universal scaling relations 
[10-12]. On the other hand, they help to study the observable properties of DM in halos of 
galaxies and clusters of galaxies [6, 13, 14] and are used in the recent numerical simulations 
[15]. These profiles also represent a very useful tool for modeling the mass distribution in 
halos [16] as well as for analytical analysis of properties of DM halos [17]. 

Galaxies, groups and clusters of galaxies are the inhomogeneities we observe inside the 
cell of uniformity which is of the order of 190 Mpc [18]. Deep inside this cell, the Universe 
is highly inhomogeneous and is well described by the discrete cosmology approach [18, 19]. 
The structure and dynamics of inhomogeneities are influenced by two main factors. They are 
the gravitational attraction between the constituents of these objects and the cosmological 
expansion of the Universe. Obviously, there is a distance from the center of mass of an 
inhomogeneity at which the cosmological expansion begins to prevail over the gravitational 
attraction. Moreover, because the expansion in the late Universe is accelerating, both of these 
phenomena act against each other. This means that the corresponding forces are directed 
oppositely. This effect was observed experimentally [20-23] and the corresponding distance 
was called the radius of the zero velocity surface. In [24] this distance has received the name 
of the radius of zero gravity. In our opinion, this distance is more properly termed as the 
radius of the zero acceleration [18, 19] because gravity does not disappear at this distance 
but the acceleration of a test body (a dwarf galaxy) is equal (or approximately equal) to 
zero. This distance for the Local Group was estimated in [25]. 

It is natural to suppose that the edge of the DM halo in clusters of galaxies corresponds 
to the surface of the zero acceleration. This idea was discussed in [17]. Here, the authors 
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considered the background metric in the form of the Schwarzschild-de Sitter (SdS) solution. 
Such approach disregards the presence of matter in the Universe. However, in the ACDM 
model matter contributes 31% into the total mass-energy balance. Therefore, in the present 
paper we include matter into account considering the Friedmann-Robertson-Walker (FRW) 
metric as the background one. 

In section 2, we start from the investigation of the gravitational potential produced 
by a spherically distributed DM halo. This potential satisfies the Poisson equation [19]. 
We supplement this equation with the proper boundary conditions at the surface of zero 
acceleration. Then, with the help of the virial theorem, we define the mean velocity dispersion 
of the galaxies which we consider as test bodies in the cluster DM halo. Here we use the 
observation [26] that distribution of galaxies in rich clusters is abruptly bounded at some 
distances much smaller than the size of a cluster. Thus, we introduce an effective radius 
Reff and suppose that all galaxies are concentrated inside a sphere of this radius. In section 
3, we apply the derived formulae to the case of the Navarro-Frenk-White (NFW) DM profile 
for the Coma cluster. Here, first, we demonstrate that Reff is close to the observable value 
if we use the NFW profile parameters for the Coma cluster found from the observations 
and simulations. Second, the comparison of our toy model with the observable data and 
simulations enables us to get preferable NFW prohle parameters for the Coma cluster. The 
main results are briefly summarized in concluding section 4. 


2 Zero acceleration sphere and mean velocity dispersion 


In the framework of discrete cosmology developed in the recent series of papers [18, 19, 
25] the acceleration of any test body in the gravitational field of a spherically symmetric 
gravitationally bound system is given by the formula (see, e.g., section 4 in [19]) 


^ dR' 


( 2 . 1 ) 


where R is the distance from the center of mass of the system to a test body, a is the scale 
factor of the Universe, the overdots denote the derivatives with respect to the synchronous 
time, and the gravitational potential satisfies the Poisson equation 


A4> 


1 cP 

Rd^ 


(R4>) = AnGNPph , 


( 2 . 2 ) 


where Gj\f is the Newtonian gravitational constant and pph is the physical (not comoving!) 
rest mass density of the bound system. It is worth noting that Eq. (2.1) was obtained in 
[19] for the Universe with hyperbolic spatial topology JC = —1. However, deep inside the cell 
of uniformity (with the size of the order of 190 Mpc [18]) the Universe can be considered 
with very high accuracy as spatially flat. Therefore, we can drop the contribution of JC to 
(2.2) and consider the Laplace operator A in cartesian (physical) coordinates. It should be 
mentioned that in the Newtonian limit the similar equations (of the form (2.1) and (2.2)) 
were obtained in [27] and are used for the N-body simulations [28]. 

The term {a/a)R in (2.1) arises due to the global cosmological expansion of the Universe 
while the term with 4*(R) takes into account the gravitational attraction between a test body 
and the gravitationally bound system. According to the recent observations, our present 
Universe undergoes the stage of the accelerated expansion. For the ACDM model, this stage 
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started at the redshift 2 : ~ 0.76 [29]. It means that from this time a > 0. Therefore, Eq. (2.1) 
demonstrates that (starting from this time) there is a distance Rh from the center of mass of 
the system at which the gravitationai attraction and the cosmoiogicai expansion compensate 
each other: 


R 


Rh 


= 0 


Rh = 


d<^/dR 


a/a 


Rh 


(2.3) 


Obviousiy, the Hubbie flows begin to form at distances greater than Rh- We can aiso rewrite 
Eq. (2.1) in the form 


9$ 


dR ’ 


(2.4) 


where 

m) = 

Then, the zero acceieration condition is 


-l-R^ + HR)- 

2 a 


(2.5) 


(9$ 


= 0 . 


Rh 


( 2 . 6 ) 


In the present paper, we consider gravitationally bound systems consisting of a rather 
large number of inhomogeneities/galaxies which have a common (for a whole system) DM 
haloes. DM halos make a dominant contribution to the total mass of the systems and 
galaxies are considered as test bodies. Such situation takes place in rich clusters of galaxies, 
e.g., for the Virgo cluster which consists of about 2000 galaxies [23, 30]. It is well known 
that the shape of this cluster is rather far from spherically symmetric [31, 32]. The cluster 
is an aggregate of at least three separate subclumps. In general, there is no problem to 
define a surface (or an approximate surface) of zero acceleration in the case of a number of 
gravitating sources which make the dominant contributions to the gravitational acceleration 
[25]. However, the DM profiles for galaxies and clusters of galaxies are usually modeled by 
the spherically symmetric distributions of DM. The Coma cluster (Abell 1656) has more 
spherical distribution^. Therefore, in what follows we consider the Coma cluster which, in a 
rough approximation, has an overall DM halo of the spherically symmetric form. 

It is natural to suppose that the edge/cutoff of a DM halo coincides with the surface 
of zero acceleration. DM particles are confined by the gravitational interaction inside the 
domain bounded by this surface and are “blown” by the cosmological expansion outside it. 
According to the observations, DM makes the main contribution into the mass of galaxies, 
groups and clusters of galaxies. Then, the rest mass density pph in Eq. (2.2) is mainly defined 
by DM. For simplicity neglecting the contribution of the visible matter as compared with 
that of DM, for the halo mass we obtain 


Rh 

M = Att J ppY,[R)R?dR. (2.7) 

0 

It is clear that the gravitational potential should satisfy the following boundary condition: 
—)■ —GnM/Rh for R —)■ Rh- Therefore, the boundary conditions for the potential <I> are: 


Im 


= 0 , 

R=Rh 


HRh) 


^Although here we also observe subhalos [33]. 
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For the given above boundary conditions, the radius of the zero acceleration sphere is 


Rh 


GnM 

aja 


-,1/3 


^{Rh) 


3 GnM 
2 Rh 


(2.9) 


Obviously, the value of Rh depends on time. On the other hand, Eq. (2.9) results in a useful 
relation 

-R% = GnM. ( 2 . 10 ) 

a 

The ratio d/a can be obtained from the second Friedmann equation. In the case of the 
ACDM model we have: 


d 

o 


- 4 
npc 
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Ac2 


= m 



( 2 . 11 ) 


where k = SttGn/c^ (with c being the speed of light) and we took into account the late time 
acceleration of the Universe expansion. In accordance with the conventional ACDM model, 
the Universe is supposed to be filled with the nonrelativistic matter (dust) characterized by 
the average rest mass density p (being constant in the comoving reference frame) and the 
dark energy represented by the cosmological constant A. We also introduce the standard 
density parameters 


= 


- 4 
KpC 


3H^a'Q 


3 ’ 


Da = 


Ac2 


( 2 . 12 ) 


where oq and Hq are the current values of the scale factor and the Hubble parameter H (t) = 
d/a, respectively. Therefore, from (2.9) we get 


Rh 


(GwM)V3 


Rn ( + 


1/3 


(2.13) 


If instead of the discrete cosmology approach the Schwarzschild-de Sitter one is applied 
(see, e.g., [17]), then in the above and following formulae we should simply put Dm = 0, 
thus disregarding the contribution of matter. However, according to the resent observations, 
matter contributes approximately 30% into the total balance (dark energy/A-term gives 
another 70%). Hence, dropping the contribution of matter corresponds to a decrease in 
accuracy of 30 percent. At the present, we are entering the era of precision cosmology. 
Therefore, future observations can reveal these discrepancies. 

The mean velocity dispersion of galaxies a in clusters is one of their important observable 
parameters [34]. With respect to the Coma cluster, the line-of-sight velocity dispersion is 
954 ± 50kms“^ [26] (the authors of [34] obtained a bit different value of cr ~ 1008 kms“^). 
We can estimate cr with the help of the virial theorem. Eq. (2.4) can be written in the form 


V - - — — 

~ ’ 

Multiplying both sides of Eq. (2.14) by R, one gets: 


|(VR)^U^-||r. 


(2.14) 


(2.15) 
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Because V and R are finite functions, averaging this equation over times which are much 
longer than characteristic periods of the system, we can eliminate the left hand side of (2.15). 
Next, we average over the volume of the halo and obtain 


— (9d> 

= —B 

dR 


(2.16) 


The mean velocity dispersion a = . If the number density of galaxies over the entire 

volume of the halo (i.e. for R £ [0, Rh]) is constant, then 

'•Rh 


- Att 

= — 
Vh 


dR 


R^dR, 


Vh = -7tR% . 


(2.17) 


However, observations indicate that distribution of galaxies is abruptly bounded at some 
distance which is considerably less that the size of the cluster. The most galaxies are 
inside this distance from the center of mass of the Coma cluster [26]. Therefore, the averaging 
over the entire volume of the halo gives the wrong result. Unfortunately, we do not know the 
real 3D picture for vectors of velocities V of all galaxies in the Coma cluster. Hence, to get 
the mean velocity dispersion we need to introduce some theoretical model. Let us suppose 
a toy model where all galaxies are concentrated inside a sphere of the radius Reff with the 
constant number density. There are no galaxies outsider this sphere. For such picture. 


— Att 
1/2 = 




eff Jo 


Reff 

^R^dR. 

dR 


^eff = -TtR: 


eff 


(2.18) 


Now we want to take a definite DM halo profile and, first, calculate the gravitational 
potential <I>, second, obtain the mean velocity dispersion from the formula (2.18) as a function 
of Reff- As we know the observable value of a as well as the experimental values for the halo 
profile parameters, we can estimate unknown Reff and compare it with the observable R 2 t- 
If these two values are close to each other, it will demonstrate that our toy model determines 
more or less successfully the characteristic radius within which most of the galaxies are 
located. On the other hand, if Reff is close to the observable value, we can use this model to 
determine parameters of the considered DM halo prohle. In the next section, we demonstrate 
our approach on the example of the NFW profile. 


3 Navarro-Prenk-White profile in the Coma cluster. Theory vs observa¬ 
tions and simulations 


At the present, the NFW profile [35] is one of the most commonly used profiles to describe 
DM distribution in galaxies and clusters of galaxies. It reads 


/5ph(7?) 


R 

R's 


4/9s 



(3.1) 


where ps = Pph{Rs) and the scale radius Rg are free parameters of the profile. Since DM 
gives the main contribution to the total mass of the system, the mass M* within a sphere of 
the radius i?* is 


M* = M{R^) 


WirpsRg 


In 



R^ 

Rg + R* 


(3.2) 
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Then, the total mass of the cluster is M = M{Rh)- Very often DM profiles are characterized 
by the mass M 200 = -Tf(i? 20 o) where the radius -R 200 defines a sphere where the average DM 
density p 200 equals 200 pcrit with 

Peril = 3Hq/{8ttGn) ~ 1.88 X 10“^®kgm“^ (3.3) 


being the critical density of the Universe today. According to the most recent Planck results 
[5], Qm ~ 0.3 and Da ~ 0.7. Besides, Hq ^ 67.8 kms“^ Mpc“^, therefore, h « 0.678. Hence, 
by the definition 

4 0 _ 

M 2 QQ = ■^T^R 2 mP '200 ) P 200 = 200pcrit . (3.4) 

It follows from Eq. (2.13) that the total mass 


87r 


1 . 


T7 — PcrilR-H ( T Da 1 — 




where 


.. _ Rh _ R 200 

= '““rT 


(3.6) 


The parameter c is usually called the concentration parameter. On the other hand, from Eq. 
(3.2) we have 

1"(1 + - irfe 


M = M 200 - 1 n , N C 
ln(l + c) - ^ 

Equating (3.5) and (3.7), we get the following useful relation: 


1 

Too 


~7:Dm + Da 


ln(l + c) — 


1 


ln(l + Xh) - 


Ah 


1 + Ah 


(3.7) 


(3.8) 


l + c| 

The ranges of parameters for the Coma cluster halo according to the observations are^ [36]: 

= 1.77/i-^Mpc = 2.61 Mpe, = 2.20 Mpe = 3.24 Mpe , 

fo, (3.9) 


= 1.32/i"^ X lO^^M 


o J 

= 2 , 


= 2.53/i"^ X lO^^Mfr 


= 17. 


Let us turn now to the gravitational potential and the velocity dispersion. To get the 
expression for the gravitational potential $ defined in (2.5), we should first solve the Poisson 
equation (2.2) and, then, use the boundary conditions (2.8). As a result, we obtain: 


HR) = -- 


R^ 9 Rh\ 

— - Rjj + ) + WttGnP. 
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Rg + R G]\jM 
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Rh {2Rl 
1 + R/Rg 
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+ In 


R In (1 + Ah) ~ i+Ag 1 + Ah 1 + Ah 


(3.10) 


^It is worth noting that if we calculate formula (3.4) with from (3.9), then 

we get a little bit different values: = 1.29/i“^ x 10 ^®Mq and = 2A7h~^ x IO^^Mq. 
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where for aja we used Eq. (2.10) and ps is expressed with the help of Eq. (3.2) for = Rji- 


Now, we can calculate the velocity dispersion with the help of Eq. (2.18). Instead of 
it makes sense to consider the dimensionless quantity a\ 



a = 


E2. 


-|l/2 


GnM 


^effj 
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5 \ Xh 


3 

2 + 
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X 


3 Xh yi y y-L,/, x s ^e// 

21 +\h ^2A^“4“ 2^ ^ eff) + ^l^ 
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where 




-^e// 
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On the other hand, from the dehnition of a we get 


CJ = 


'e2. Rh 


GnM 


1/2 


E2 


1/2 


SVT Gn X‘jjR 2 QQP crit 




-n 1/2 


(3.11) 


(3.12) 


(3.13) 


where Rh is expressed as Rh = R20oXh/c and the total mass M is given by Eq. (3.5) with 
M 200 from (3.4). 

According to the observations [26], the line-of-sight velocity dispersion for the Coma 
cluster is 954 ± 50kms“i. To get the 3D velocity dispersion we should multiply this value 

Tl 1/2 


by \/3. Therefore, the minimal and maximal 3D velocity dispersions are: (E(™”'))2 


1566 km s ^ and 


.nl/2 


maximal values of a are: 

^(rain) _ 


(E(™-“*))2 = 1739kms 1 , respectively. Therefore, the minimal and 


(T' 


{y+ni+Y 


1/2 




-11 1/2 


(3.14) 


cr 


{rnax) _ 


(l/(maa;))2 


■ 11/2 




-1 


1/2 


(3.15) 


where 72200 "^’”^°^*^ should be taken from (3.9). 

According to Eq. (3.9), for the Coma cluster the concentration parameter c € [2,17]. 
Our task is to determine the values of Reff as a function of c from this region. To perform it, 
we can start from c = 2 and then define Reff for each value of c with the step, e.g., Ac = 1 up 
to c = 17. It is convenient to present the obtained values of Reff in the corresponding table 
Reff vs c. To demonstrate this procedure, let us consider as an example the case c = 4. Eor 
a given value of c, the parameter Xh can be found from Eq. (3.8). In the considered example 
this is Xh = 33.28. For this value of Xh, the graph of a (see Eq. (3.11)) as a function of Ag// 
is drawn in Fig. 1. From Eqs. (3.14) and (3.15), we can also get for c = 4: = 1.15 

(lower horizontal line in Fig.l) and = 1,60 (upper horizontal line in Fig.l). It can be 

easily seen from this picture that for each value of have two roots. They are 

Xeff,i = 0.39 and Ae //,2 = 21.31 for Ae //,3 = 1.58 and Ae //,4 = 5.68 for Since 
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^effA 

4.17 

3.70 

3.71 

3.72 

3.74 

3.75 

3.76 

3.77 

c 

2 

3 

4 

5 

6 

7 

8 

9 

^effA 

3.78 

3.78 

3.79 

3.80 

3.80 

3.80 

3.81 

3.81 

c 

10 

11 

12 

13 

14 

15 

16 

17 


Table 1. The values of the root for different values of the concentration parameter c. 


^{rmn) ^^{max)^ obtained with the help of corresponding values 

for Rh = R 200 ^h/c and Rg = R 200 /C are = 26.99 Mpc and = 0 . 8 I Mpc 

(^^in) _ 21/12 Mpc and = o.65 Mpc). Therefore, from the relation Reff = XeffRs 

we obtain four marginal values: = 0.32 Mpc and R^^/^ = 17.28 Mpc (for 

and Ae//i, 2 ) correspondingly) and = 1-03 Mpc and = 3.71 Mpc (for 

and Ae// 3 , 4 , correspondingly). 

These calculations demonstrate that only is close to the observable value R^t ~ 3 

Mpc [26]. Calculations for other values of the concentration parameter c from the considered 
region lead to the same conclusion. Hence, only the marginal value corresponds to 

the observations. We summarize the results of our calculations in Table 1. 


a 



Figure 1. This plot shows the behavior of the dimensionless velocity dispersion ct as a function of 
Ae// for the concentration parameter c = 4 (curved line). Upper and lower horizontal lines correspond 
to the maximal and minimal values of a obtained with the help of the observable data. 


This table shows that the closest with respect to the observations results take place for 
the values c = 3 4- 4 where has a minimum. 

It is worth noting that this concentration parameter is very close to the one obtained 
from the fitting of the NFW model. This fit is based on a series of DM cosmological simula¬ 
tions within the ACDM model [15, 37]. For example, if we take = 1.29h~^ x 

(see footnote 1), then we get from Eq. (24) in [15] or equivalently Eq. ( 8 ) in [37] that 
c « 3.89. 

Therefore, the comparison of our toy model with the observable data and simula¬ 
tions leads to the following preferable NEW parameters for the Coma cluster: R 200 ~ 
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1.77h~^ Mpc = 2.61 Mpc and c = 3 h- 4. For these values of c and M 200 = 1.29/i“^ x IO^^Mq, 
the total mass of the Coma cluster in the case of the NFW profile is M = (4.08^-4.74) h~^ x 
IO^^Mq (see Eq. (3.7)). We also found that in this cluster the most of galaxies are concen¬ 
trated inside a sphere of the effective radius Reff ~ 3.7 Mpc and the line-of-sight velocity 
dispersion is 1004 km s“^. This value is very close to 1008 km s“^ obtained from the observa¬ 
tions for the Coma cluster (Abell 1656) in [34]. 

4 Conclusion 

In the present paper we have investigated the distribution of visible and dark matter in the 
Coma cluster. In the late Universe, such clusters are gravitationally bound systems and can 
be considered as isolated objects. To study them, we can apply the discrete cosmology (me¬ 
chanical) approach [18, 19]. In this approach, the motion of a test body in the vicinity of an 
inhomogeneity is defined by two factors: the gravitating attraction of the inhomogeneity and 
the accelerated cosmological expansion of the Universe. These two competitive mechanisms 
define a region where the acceleration of a test body is equal to zero; the zero acceleration 
surface. In the case of spherical distribution of matter this surface is a sphere of the radius 
Rh- Inside this sphere, the gravitational attraction plays the main role while at the larger 
distances the cosmological expansion prevails. This is the region of the Hubble flows forma¬ 
tion. It is well known that the DM halo gives the main contribution to the total mass of a 
cluster. Hence, the gravitational potential of a cluster is mainly defined by the DM distri¬ 
bution in it. There are a number of halo profiles proposed to describe the DM distribution 
in galaxies and clusters of galaxies. It is natural to define the edge of the DM halo at the 
radius of zero acceleration. Because the halo of DM dominates in clusters of galaxies, the 
visible matter in the form of galaxies is considered as test bodies in the gravitational held of 
the halo. 

In our paper, we have investigated the Poisson equation for the gravitational potential 
of a spherically symmetric halo prohle. We have dehned the proper boundary conditions 
at Rh- Then, with the help of the virial theorem, we have got the formula for the mean 
velocity dispersion for the visible matter (galaxies) inside this halo. Observations show [26] 
that galaxies inside clusters are abruptly bounded at some distance which is considerably 
less that the size of the cluster. The most galaxies are inside this distance from the center of 
mass of the Coma cluster [26]. Since we do not know the real 3D velocities of galaxies, we 
have supposed a toy model where all galaxies are concentrated inside a sphere of an effective 
radius Reff with the constant number density. This has enabled us to get the expressions 
for the mean velocity dispersion for an arbitrary spherically symmetric halo profile. 

To perform specific calculations, we need to suggest a definite form for the DM density 
profile. The NFW profile [35] is one of the most commonly used. We have used this profile in 
the case of the Coma cluster. First, we have got the formula for the gravitational potential 
in this profile. Then, we have obtained the mean velocity dispersion as a function of Reff- 
We have shown that the calculated value of Reff can be rather close to the observable 
R 2 t within the observation accuracy of the NFW parameters. Moreover, the comparison of 
our toy model with the observable data and simulations leads to the following preferable 
NFW parameters for the Coma cluster: R200 ~ 1-77 x h~^ Mpc = 2.61 Mpc, c = 3 4 and 

-^^200 = l- 29 h~^ X IO^^Mq. We have also found that in this cluster the most of galaxies are 
concentrated inside a sphere of the effective radius Reff ~ 3.7 Mpc, and the line-of-sight 
velocity dispersion is 1004kms“^. 
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